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ABSTRACT 



Aims. The line profile variability and photometric variability of the 09.5 Vp star HD 93521 are examined in order to establish the 
properties of the non-radial pulsations in this star. 

Methods. Fourier techniques are used to characterize the modulations of the He i AA 5876, 6678 and H a lines in several spectroscopic 
time series and to search for variations in a photometric time series. 

Results. Our spectroscopic data confirm the existence of two periods of 1.75 and 2.89 hr. The line profiles, especially those affected 
by emission wings, exhibit also modulations on longer time scales, but these are epoch-dependent and change from line to line. Unlike 
previous claims, we find no unambiguous signature of the rotational period in our data, nor of a third pulsation period (corresponding 
to a frequency of 2.66 d -1 ). 

Conclusions. HD 93521 very likely exhibits non-radial pulsations with periods of 1.75 and 2.89 hr with / 8 ± 1 and / - 4 ± 1 
respectively. No significant signal is found in the first harmonics of these two periods. The 2.89 hr mode is seen at all epochs and in 
all lines investigated, while the visibility of the 1.75 hr mode is clearly epoch dependent. Whilst light variations are detected, their 
connection to these periodicities is not straightforward. 
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1. Introduction 

Time resolved spectroscopy of O-type stars has shown that 
absorption line profile variability at the level of a few per cent 
is a common feature (e.g. Fullerton et al. 1996 ). Various mech- 
anisms, including magnetic fields, stochastic structures at the 
base of the wind and non-radial pulsations have been proposed 
to explain this variability. Despite their rather subtle signature, 
the diagnostic potential of these phenomena is considerable. 
Especially in the case of non-radial pulsations, the emerging 
discipline of asteroseismology offers the possibility to gain 
insight into the interiors of early-type stars. However, to charac- 
terize the nature of the phenomenon requires a rather long and 
well sampled time series of spectra with a high resolution and a 
high S/N (see e.g. Aerts et al. l2004l for the case of the /JCephei 
variable vEri). Up to now, such detailed studies have therefore 
been restricted to a few, rather bright and well-known O-stars 
such as ( Pup (Baade et al. [19911 ) and ( Oph (Kambe et al. [T997] l. 

In this context, the high Galactic latitude O-star HD 93521 
(hi = 183.14°, bu = 62.15°) is an extremely interesting target. 
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Based on optical spectroscopy, Fullerton et al. ( 119911 and 
Howarth & Reid ( fF?93j l found HD 93521 to exhibit bumps at 
the 1 % level moving from blue to red across the profiles of the 
He i lines, whilst no variability was detected in the Hen lines. 
Fullerton et al. ( 119911 ) as well as Howarth & Reid d 19931 ) accord- 
ingly interpreted these features as the signature of non-radial 
pulsations with a period of 1 .8 hours. The He i - He n dichotomy 
was interpreted as arising from the substantial gravity darkening 
that favours He i line formation near the cooler equatorial regions 
where the pulsational amplitude attains a maximum (Townsend 
119971 ). The existence of non-radial pulsations was subsequently 
confirmed by Howarth et al. ( 1998 ) using IUE spectra and three 
different periods were identified. Recently, Rzaev & Panchuk 
(2006) reported on the existence of slightly different variability 
patterns between the strong and weak Hei lines. However, since 
the Rzaev & Panchuk (2006) data set covered only 2.7 hours, 
it yields hardly any constraint on the properties of the pulsations. 

HD 93521 has one of the largest rotational velocities known 
among O-stars (34 1 km s~ 1 , Penny 1 996 ; 432 km s~ 1 , Howarth et 
al. 119971 390kms~ 1 , see below in this paper). The stellar wind 
has an apparently low terminal velocity and is likely heavily dis- 
torted by rotation (e.g. Bjorkman et al. 1994). In the optical, the 
wind produces emission features in the wings of the Ha line, al- 
though this line has so far never been reported as a pure emission 
feature. 

While the optical spectrum of HD 93521 leads to an 09.5 V 
classification, the nature of this star has been subject to debate 
over many years. In fact, assuming a typical absolute magnitude 
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Table 1. Summary of our observing runs (see text for the meaning of the different columns). The green and red wavelength ranges 
stand respectively for 5500 - 5920 and 6530 - 6710 A, whilst SPM stands for the echelle spectra taken at San Pedro Martir. The last 
column yields the mean signal-to-noise ratio of the spectra gathered during the corresponding campaign. The last line refers to the 
photometric monitoring campaign that took place in coordination with the April 2005 spectroscopic campaign. 



Epoch 


Domain 


AT (days) 


N 


At (days) 


AVna, (d- 1 ) 


Vmax (d~') 


S/N 


February 1997 


red 


4.155 


35 


2.49 x 10~ 2 


0.241 


20.1 


165 


April 2005 


red 


3.809 


55 


1.72 x 10~ 2 


0.263 


29.1 


280 


April 2005 


SPM 


5.096 


90 


0.84 x lO- 2 


0.196 


59.6 


250 


February 2006 


green 


6.173 


75 


1.60 x 10 2 


0.162 


31.3 


650 


April 2006 


green 


0.180 


13 


1.49 x 10~ 2 


5.55 


33.5 


580 


April 2007 


green 


6.191 


56 


1.46 x 10~ 2 


0.162 


34.2 


550 


April 2005 


U BB { B 2 VViG 


28.221 


378 


0.31 x IQ 2 


0.035 


161.3 






Fig. 1. Average spectrum (top) and temporal variance spectrum (TVS, bottom) computed from the green, red and SPM spectra 
taken in February + April 2006 (left panel), April 2005 (middle) and April 2005 (right) respectively. The dotted lines yield the 99% 
significance level for the variability evaluated following the approach of Fullerton et al. (1 19961 1. 



for a Population I 09.5 V star, HD 93 521 i s located about 
1.4kpc above the Galactic plane (Irvine 1989), far away from 
any site of recent massive star formation. While there is still 
some uncertainty concerning the motion of the star (towards or 
away from the Galactic plane, see Gies ll987l . it seems unlikely 
that HD 93521 could have formed in the plane and subsequently 
moved to its current position. Furthermore, based on the rather 
low stellar wind terminal velocity of HD 93521 and assuming 
that Voo = 3 x v esc , Ebbets & Savage (|1982| l concluded that 
this star was likely a low-mass evolved Population II object. 
However, Irvine (1989) showed that the relation between Vco 
and v esc does not hold for late O-type stars such as HD 93521 
and he proposed that the star is in fact a normal main sequence 
star that has formed in the Galactic halo (another example of a 
massive star born in the halo can be found in Heber et al. 1995). 
Lennon et al. (1991) measured the equivalent widths (EWs) of 
several metal lines in the spectrum of HD 93521. Though these 
lines are washed out by rotational broadening and the EWs are 
affected by large uncertainties, their strengths are inconsistent 
with Population II metal abundances. Finally, Massa (1995) 
showed that the other peculiar features of this star (its unusually 
low UV continuum, its abnormally strong wind lines in the UV, 
its low excitation photospheric lines) can all be accounted for 
by the effect of gravity darkening in a 'normal' Population I star 
rotating at 90% of its breakup velocity and seen nearly equator 
on. 



In this paper, we present the results of a spectroscopic and 
photometric monitoring campaign of HD 93521. The aim of this 
project was to check the long-term stability of the periodicities 
identified in previous investigations. 

2. Observations 

2.1. Spectroscopy 

We observed HD 93521 during five observing campaigns (in 
February 1997, April 2005, February and April 2006 as well as 
April 2007) with the Aurelie spectrograph at the 1.52 m tele- 
scope of the Observatoire de Haute Provence (OHP). The 1997 
and 2005 data were taken with a 1200 lines/mm grating blazed at 
5000 A and covered the wavelength domain from 6530 to 6710A 
with a resolving power of R 21000. The 2006 and 2007 data 
were obtained with a 600 lines/mm grating also blazed at 5000 A 
and covering this time the region between 5500 and 5920 A with 
R - 10000. In 1997, the detector was a Thomson TH7832 linear 
array with 2048 pixels of 13 /im each, whilst later observations 
used an EEV 42-20 CCD with 2048 x 1024 pixels. All th e OHP 
data were reduced in the standard way (see Rauw et al. 2003) 
using the MIDAS software provided by ESO. To correct the ob- 
served spectra to first order for telluric (mainly water vapor) ab- 
sorption lines, a template of the telluric spectrum was built from 
observations of a bright star at very different airmasses (usually 
around 1.2 and 3.0). The stars used for this purpose were vOri 
(B3 V) and HD 100889 (B9.5 Vn) in February 1997, HD 177724 



G. Rauw et al.: Spectroscopic and photometric variability ofthe 09.5 Vp star HD 93521 



3 



(AO Vn) in April 2005 and Regulus (a Leo, B7 V) in February 
2006. Finally, the spectra were normalized using properly cho- 
sen continuum windows. Special care was taken to ensure that 
all the data were normalized in a homogeneous way to allow a 
self-consistent search for variability. 

In April 2005, we gathered a series of 90 spectra with the 
echelle spectrograph mounted on the 2.1m telescope at the 
Observatorio Astronomico Nacional of San Pedro Martir (SPM) 
in Mexico. The instrument covers the spectral domain between 
about 3800 and 6800 A with a resolving power of 18000 at 
5000 A. The detector was a SITe CCD with 1024 x 1024 pix- 
els of 24 /mi 2 . The data were reduced using the echelle context 
of the MIDAS software and specific orders covering several im- 
portant lines (He i AA 44-7 1, 5876, HB and Ha) were normalized 
using a set of selected continuum windows. Unfortunately, the 
He i A 6678 line could not be studied because it fell too close to 
the edge of an order. 

A summary of the main characteristics of our various data 
sets is provided in TableQ] For each run, AT gives the time 
elapsed between our first and our last observation, while N is 
the total number of observations. At is the average time inter- 
val between two consecutive exposures during the same night. 
In light of the Fourier analysis in Section [3~2l we further list the 
typical FWHM of a peak in the power spectrum A v nat = 1/AT, 
as well as v max = 1/(2 Af) which provides a rough indication 
of the highest frequencies that could be sampled with our time 
series. 

During the April 2006 campaign, HD 93521 was also ob- 
served in the near-infrared with the Aurelie spectrograph. The 
data were obtained with a 300 lines/mm grating blazed at 6000 A 
and covering the region between 8055 and 8965 A with R 
6500. Whilst the individual spectra in this wavelength domain 
do not have a sufficient quality to perform a meaningful vari- 
ability analysis, we present the mean spectrum over the useful 
band from 8420 to 8900 A in Fig. [2] The spectrum is largely 
dominated by broad and shallow hydrogen absorption lines of 
the Paschen series (from Pal 1 to at least Pal6) with some mod- 
est contributions from Hei and possibly Hen, Cm and Nm 
lines. The RVS instrument onboard ESA's forthcoming astrome- 
try mission GAIA will cover a sub-domain (from 8470 - 8740 A) 
of the band shown in Fig.|2]with a resolving power twice as large 
as that of our data to measure among other things the radial and 
rotational velocities. Our spectrum indicates that the only promi- 
nent features in the RVS spectra of rapidly rotating late O-type 
dwarfs are the Paschen lines. 




8500 8600 8700 8800 8900 

Wavelength (A) 



Fig. 2. Mean spectrum of HD 93521 over the 8420 - 8900 A 
wavelength domain. The ions responsible for the lines (either 
stellar, or telluric in the case of the OH emissions) are indicated 
above the spectrum, whilst the diffuse interstellar bands (DIBs) 
are indicated below. 



2.2. Photometry 

Time resolved photometry of HD 93521 through the 
(U B B\ B2VV\G)c Geneva filter system (see e.g. Rufener 
119811 Bessell [2005) was obtained with the Mercator telescope 
in April 2005. Mercator is a 1.2 m semi-robotic telescope, 
situated at Roque de los Muchachos (La Palma Island, Spain). 
The telescope is equipped with the P 7 photometer, which is a 
two channel photometer for quasi simultaneous 7 band mea- 
surements. The first channel (A) is centered on the star while 
the second channel (B) is centered on the sky. The position 
angle of the sky can be changed by turning the derotator while 
changing the distance between both channels needs manual 
interaction. The filter wheel turns at 4 Hz and a chopper directs 
both channels alternatively to the photomultiplier. As such, 
the photomultiplier measures both beams A and B through the 
seven filters four times each second. 

The strategy for performing the observations is oriented to- 
wards obtaining high-precision photometry. In order to achieve 
this, we measure stars within a range of 0.1 in air-mass F z 
for nights of good atmospheric conditions. It is advantageous 
to measure a variable star with a period of the order of hours 
frequently each night, so we have opted for the interval F z e 
[1.05; 1.15] to measure HD 93521. For the reduction process of 
such a type of observing night, typically 2-3 standard stars of 
different colour are observed each hour. The determination of 
the extinction coefficients was done according to the method out- 
lined in Burki et al. (1 19951 1 and results in measurements with a 
typical accuracy of a few mmag for the colours of stars with vi- 
sual magnitude between 5 and 10. A total of 378 data points were 
gathered between 29 March and 27 April (UT dates). In our anal- 
ysis we have restricted ourselves to those 315 data points that 
have weights of 3 or 4 (for good or exceptional quality measure- 
ments respectively, see Rufener [l9811 l. Nights with these quality 
weights are characterized by standard deviations on the measure- 
ments of standard stars of less than 0.005 mag. 



3. Results 

3.1. Spectroscopy 

The average spectra of the star over the wavelength domains in- 
vestigated here are shown in Fig.Q] Besides a few diffuse inter- 
stellar bands (DIBs), identified following the catalogue of inter- 
stellar features from Herbig (1995), we note the broad absorption 
lines as well as the emission components of Ha and He i A 5876. 
Note that the He i lines appear relatively strong compared to the 
Ha line. This can be understood in view of the helium overabun- 
dance (y - 0.18 ± 0.03) reported by Lennon et al. ( 119911 ) and 
Howarth & Smith d200Tb . 

In addition to the Om i5592, Civ ^5801, 5812 and Hei 
A 5876 absorptions that are typical for an O-type star, the spec- 
trum in the region between 5500 and 5900 A also reveals several 
unusual features, such as the N n AA 5667, 5676, 5680 blend and 
the Si m A 5740 line. These lines are strong in the spectra of early 
B supergiants, but are usually absent from O-star spectra (see 
Walborn 1980). The same remark applies to the feature near the 
red emission wing of Ha that we identify as Cn /LI 6578, 6583 
which are clearly seen in the spectra of early B-type stars. These 
particularities are reminiscent of the description of the IUE spec- 
trum of the star given by Massa ( 1995) and thus likely reflect the 
strong temperature gradient that exists at the stellar surface as a 
result of gravity darkening. According to Massa ( 119951 ), the po- 
lar region has a temperature close to that of an 07 star, whilst the 
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HD 93521: He I X 6678 



X 
2 

CD 

.SO 

go 



1 

98 
96 
94 
92 




6670 6680 

Wavelength (A) 



6690 



- 

o 
_ 

O 
J 




1 -0.5 
Log(a (A" 1 )) 



Fig. 3. Determination of the rotational velocity of HD 93521 us- 
ing the Fourier technique. The top panel illustrates the mean pro- 
file of the He i A 6678 line as observed in April 2005. The lower 
panel illustrates the resulting Fourier transform along with the 
Fourier transform G{cr) of a rotational broadening function (see 
Gray 2005 ) evaluated for v sin i = 395 km s . 



equatorial region is much cooler and experiences a lower surface 
gravity, hence producing a spectrum similar to that of an early B 
supergiant. 

We have determined the star's projected rotational velocity 
v sin / applying the Fourier method (see Gray 2005 Simon-Diaz 
& Herrero 120071 1 to the mean profiles of the Hei AA5876, 
6678 and Oin A 5592 lines. To first order, the structure of the 
Fourier transforms of these lines can be explained by projected 
equatorial rotational velocities of 345, 395 and 380km s _1 
respectively. While the results for the Hei A 6678 and Ora 
A 5592 line are in very reasonable agreement, Hei /1 5876 
yields a significantly lower value. There are several reasons 
why the three lines do not yield exactly the same result. First 
of all, the profile of Hei A 5876 is affected by emission that 
could fill-in the absorption in the wings of the line, hence 
leading to an apparently lower v sin i. In addition, although we 
have used the mean profiles evaluated from a large number of 
observations, the line profile variability (see below) leaves some 
residual structure in the profiles. These features lead in turn to a 
Fourier transform that deviates from the one expected for a pure 
rotational broadening^]. The cleanest result is obtained for the 
Hei A 6678 line (see Fig.[3}- In summary, we conclude that the 
projected equatorial rotational velocity of the star is most likely 
390 + lOkms 1 . There is a caveat here: our determination of 
v sin i likely represents a lower limit to the actual value. Indeed, 
Townsend et al. (2004) argued that the rotational velocities 
of very rapidly rotating Be-type stars might be systematically 
underestimated as a result of the effect of gravity darkening that 
would reduce the weight of the equatorial region (and hence 



of the highest rotational velocities) in the formation of the 
absorption line profiles. For B0 stars, Townsend et al. find that 
the He i A 447 1 and Mg n A 448 1 lines underestimate v sin i by 
about 10%. A similar underestimate could hold in our case. 

To quantify the variability of the spectral lines, we have 
computed the temporal variance spectrum (TVS, Fullerton et al. 
1996) for the data of each observing campaign. The TVS dia- 
grams of the spectral regions investigated here reveal no signif- 
icant variability in the metallic lines, except perhaps for some 
marginal variations in the region between 5660 and 58 10 A. The 
significant variability is restricted to the Hei A 4471, H/3, Hei 
A 5876, Ha and He i A 6678 lines. The majority of these lines (ex- 
cept perhaps the H/3 line, see Fig.[T} exhibit evidence for emis- 
sion wings. 



3.2. Fourier analysis 

For each of the variable lines, we have performed a 2-D Fourier 
analysis (see e.g. Rauw et al. 120011 [2003) where the Fourier 
power spectrum is calculated at each wavelength step across the 
line profile using the technique of Heck et al. 0985) revised by 
Gosset et al. ((2001). This algorithm is specifically designed to 
handle time series with an uneven sampling. The periodograms 
were calculated for each line and for each campaign up to a max- 
imum frequency of 30 d~' . 

For the frequencies that are identified from the peri- 
odograms, we have further applied the different techniques out- 
lined in Rauw et al. (2001 ). For each line and for each observing 
campaign, we have fitted an expression of the type 



I(A, t) = I (A) + J] AiOl) sin [2 nv it + 



(1) 



to the time series. Here I(A, t) is the line intensity at wavelength 
A and at time f, v, are the frequencies corresponding to the high- 
est peaks in the periodogram, whilst Af(A) and 4>i(A) are respec- 
tively the semiamplitudes and the phases for these frequencies as 
a function of wavelength (see Rauw et al. 12001 B The number q 
of frequencies that were simultaneously fitted was progressively 
increased from 1 up to 3. This method not only allows us to char- 
acterize the properties of a given frequency, but further enables 
us to 'prewhiten' the data by removing the modulation at this 
frequency from the time series and re-compute the periodogram 
for the 'prewhitened' data set. We start by removing the high- 
est peak frequency found in the periodogram of the original data 
set. Next, we identify the highest frequency in the periodogram 
of the new (prewhitened) time series and prewhiten then the 
data simultaneously for the two frequencies. The procedure is 
repeated for up to three different frequenciefl Indeed, three fre- 
quencies are usually sufficient to reduce the residual amplitude 
in the prewhitened Fourier power spectra to a level compatible 
with the noise level of the data. 

To evaluate the error bars on the semiamplitudes and phases 
of the various Fourier components, we have used Monte Carlo 



1 Macroturbulence also affects the shape of the Fourier transform, 
but this is mainly an issue for very slow rotators which is not the case 
of HD 93521. 



2 These techniques are especially useful if the line profile variations 
are interpreted in terms of non-radial pulsations. Indeed, in this case, 
Telting & Schrijvers ( 1997 ) and Schrijvers & Telting ( 1999 ) have shown 
that the observable phase differences between the blue and red line 
wings can be related directly to the degree / and to the absolute value of 
the azimuthal order \m\ of the pulsation mode. 

3 We caution that some of the low frequencies identified throughout 
this work likely correspond to long-term trends rather than to actual 
periods. 
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April 2005: Ho I A 6678 
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Fig. 4. Left: the top panel yields the mean Fourier power spectrum of the February 1997 Hei A 6678 data as derived from the 2- 
D Fourier spectrum averaged over the profile of the He i A 6678 line (from 6660 to 6700 A). The middle panel indicates the mean 
power spectrum prewhitened for the two dominant frequencies (0.54 and 7.30d~') simultaneously. The lower panel yields the power 
spectral window resulting from the sampling of the February 1997 time series. Right: same for the April 2005 data. The frequencies 
used for the prewhitening are v\ = 13.68 and v-i - 8.31 d _1 in this case. 



simulations assuming that the uncertainties on I(A, t) are equal 
to the noise level <x(?) of the corresponding spectrum evaluated 
in nearby continuum windows. In a subsequent step, we have 
also included the first harmonics of the three most significant fre- 
quencies (hence accounting for six frequencies in Eq.[T|i. Usually 
however, no significant power was found in the first harmonics 
of those frequencies that are likely associated with pulsations 
(see below). The amplitudes of these harmonics are found to be 
very low, much lower than those of the actual frequencies. The 
resulting uncertainties on the phase constants of the harmonics 
are very large, rendering any attempt to measure their blue to red 
phase shift hopeless. 

The evaluation of the absolute significance of a peak in a pe- 
riodogram of a time series with an uneven sampling is a contro- 
versial problem, that does not have an exact solution (see Heck 
et al. 1 19851 >. Therefore, as a first order indication of the signifi- 
cance of a peak, we have evaluated, at each wavelength step A, 
the quantity z v {A) = P v (A)/s v (A) 2 where P V (A) is the power in the 
periodogram at the frequency v under investigation and s v (A) 2 
is the variance of the time series after prewhitening of this fre- 
quency. Whilst z v (A) is not constant over the width of the various 
lines studied here, all the frequencies considered in the follow- 
ing yield rather large (> 20) values of z v (A) over a significant 
part of the line profile. Therefore, all the frequencies discussed 
below are detected with a good significance. However, this does 
not imply that all the frequencies reported below actually have a 
physical meaning. Some of them, especially the low frequencies, 
could rather result from transient features in the wind of the star. 



In the following, we shall refer to v\ and V2 as the 13.6 and 
8.3 d" 1 frequencies previously reported in the literature. 



3.2.1. He i ,16678 

The power spectrum calculated from the February 1997 time se- 
ries of this line reveals its highest peak at 0.54 d _1 (see Fig.|4j». 
Beside this low frequency peak, the periodogram exhibits a fam- 
ily of aliases associated with 7.30 d _1 . In the raw periodogram, 
the one-day alias at 6.30 d -1 is stronger than the one at 8.29 d _1 . 
However, this situation changes when the data are prewhitened 
for the 0.54 d _I modulation with the 8.29 d -1 peak becoming 
stronger than the one at 6.30 d -1 . 



Fig. 5. Grey-scale 2-D Fourier power spectrum of the He i A 6678 
line as derived from the April 2005 data. 

The lower frequency apparently produces a modulation of 
the entire line profile; there is no significant change of the phase 
as a function of wavelength. This variation is therefore unlikely 
to be related to non-radial pulsations. On the other hand, the 
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7.30CT 1 frequency is an alias of the V2 = 8.28 cT 1 frequency 
reported by Howarth et al. ( 119981 ). We have prewhitened the data 
using either the 7.30 d 1 or the 8.29 d -1 frequencies. Due to the 
strong level of the noise 'continuum' in the February 1997 peri- 
odogram, the prewhitening procedure does not result in a clean 
periodogram. Still, the results obtained when prewhitening ei- 
ther of the alias frequencies are the same within the level of the 
noise in the periodogram. Both candidate frequencies display a 
progression of the phase constant </>,■ from the blue wing to the red 
wing. This progression is not fully monotonic (at least partially 
because of the noise level of the data). Nevertheless, assuming 
v e sin i = 395 km s _1 , we find that the blue - red phase difference 
A <p amounts to about 3 n radians, irrespective of the choice of the 
alias. However, we stress that the amplitude associated with V2 
is largest over the radial velocity interval -350 to +100kms~ I : 
there is almost no significant signal outside this velocity domain 
(i.e. between +100kms~ I and +v e sin/). 




-500 500 




-500 500 




■s- 



-500 500 

Velocity (km/s) 

Fig. 6. The top panel yields the mean normalized Hei A 6678 
line profile as observed in April 2005 shown as a continuous 
line, as well as the mean profile after prewhitening (i.e. the Io(A) 
term in Eq.[T|) shown by a dashed line. The next panels (from 
top to bottom) yield the semiamplitudes (in units of the contin- 
uum level) of the v\ (dots, second panel) and V2 (crosses, third 
panel) modulations, the phase constant (f>\ of the v\ signal and 
the phase constant <p2 of the V2 component. The error bars were 
obtained from Monte Carlo simulations (see text). Phase 0.0 was 
arbitrarily set to the time of the first observation of this data set. 

The power spectrum of the April 2005 data looks completely 
different (Fig.|4]i. The periodogram is now dominated by a strong 
peak at v\ = 13.68 d 1 along with its one-day aliases. When 
the data are prewhitened for the v\ frequency, another family of 



aliases associated with V2 = 8.31 d 1 becomes dominant. Again, 
the mean power spectrum, averaged between 6660 and 6700 A, 
reveals some power also at lower frequencies (near 0.25 d _1 ), 
though the 2-D power spectrum indicates that this latter feature 
is mainly restricted to a narrow range of wavelengths in the line 
wings and mostly in the blue wing near 6671 A (see Fig. [5]). 

Although the power spectrum reveals some structure near 
26 d~', the frequency of the strongest peak in this region 
(25.8d _1 ) is quite different from 2xv, (= 27.36d _1 ) and it 
seems therefore unlikely that this feature is related to the first 
harmonic of v\ . The periods corresponding to v\ and V2 are re- 
spectively 1.75 and 2.89 hr, in excellent agreement with the pe- 
riods (Pi = 1. 76 + 0.04 hr, P 2 = 2.90 + 0.12hr) reported by 
Howarth et al. ( fT998b . 

Figure|6]illustrates the mean line profile as observed in April 
2005, as well as the semiamplitudes and phases of the two dom- 
inant frequencies (vi and vi) during this campaign. Again, we 
note the diagonal progression of the phase constant as a function 
of wavelength which is typical of non-radial pulsations. From 
Fig. [6] we infer a blue to red phase difference of A <pi -In and 
A <p2 — 3.5 n for the v\ and vo frequencies respectively, although 
for V2, the progression of the phase constant is not strictly mono- 
tonic across the line profile. 



3.2.2. Ha and H/3 

The power spectrum of the February 1997 Ha data shows the 
highest peak at 0.30 d -1 and its one-day alias at 0.70 d" 1 . A 
group of secondary peaks is found at 6.94 d _I and its aliases. 
We further detect some weaker peaks at 6.26 d _I (alias of vi) as 
well as its aliases. The power spectrum is rather well cleaned 
by prewhitening these three frequencies (see Fig.|7]l. The low 
frequency variation corresponds to a global modulation of the 
line intensity in the core (radial velocity interval from -100 to 
+ 100kms _I ) with a rather constant phase. Since the 6.26 d -1 
frequency is an alias of v%, we have also prewhitened the data 
using 0.30, 6.94 and 8.28 d -1 . For the latter frequency, we ob- 
serve again a blue to red phase difference of about 3 n (within 
the uncertainties). However, as for He i A 6678, the amplitude of 
this modulation is really significant only in the blue wing of the 
line. 

In the April 2005 OHP Ha data, the highest peak of the 
power spectrum occurs at 0.15 d -1 and its aliases (1.18 and 
0.78 d~'). Prewhitening either of these frequencies yields a set 
of secondary peaks near 1.67 and 3.09 d _1 . If these frequencies 
are removed in turn, the power spectrum of the residual data set 
displays some low amplitude peaks at 8.32 and 13.14d (note 
that while the former frequency is likely equal to v-i, the latter 
one is formally not consistent with vi). 

The April 2005 SPM echelle spectra yield a similar pic- 
ture: the bulk of the power of the Ha variations is again con- 
centrated in variations that occur on low-frequency timescales 
(1.09 and 1.32 d _1 , which are also found for the He i A 5876 line 
in the same data set) and affect mainly the wings of the line. 
Prewhitening the data for these two frequencies yields a high- 
est peak at 2.51 d" 1 as well as a series of aliases of v\ + 1 and 
V2 - 1. For HyS, the SPM data reveal the strongest peaks at 0.01 
and 0.44 d" 1 . Both, V2 and v\ are clearly detected though and 
the blue to red phase differences are about 3 n and 7 n, though it 
should be stressed here that the evolution of the phase constants 
is not monotonic in the case of this line. 
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Fig. 7. Same as Fig.[4]but for the He line: February 1997 (left) and OHP April 2005 (right). The power spectra are averaged over 
the wavelength domain from 6540 to 6580 A. The frequencies used for prewhitening are 0.30, 6.94 and 6.26 d -1 for the February 
1997dataand 1.18, 3.09 and 1.67c! 1 forthe April 2005 campaign. 



3.2.3. He izl/15876, 4471 

The power spectrum of the OHP 2006 data is dominated by low 
frequencies. When all the OHP data taken in February and April 
2006 are analysed together, the highest peak is found at 0.02 d _1 
(along with its aliases, Fig.[8}. Prewhitening the time series for 
this frequency and its first harmonics, yields however some arte- 
facts in the mean line profile (Iq(A) in Eq.[TJ. Hence, we rather 
used the 1-day alias 1.10 d. When we prewhiten the data for 
the latter frequency, the highest peaks in the new periodogram 
are found at 12.63 and 8.26 d -1 (= V2). The former frequency is 
likely an alias of v\. Prewhitening the frequencies 1.10, 8.26 and 
13.60d~', we find that the blue to red phase difference amounts 
to 6 n for v\ (for 12.63 d _1 , this becomes =s 7 n). For V2 we find 
a phase difference of about 3.5 n with the caveat that the phase 
progression is apparently not monotonic for this frequency (see 
Fig. [8]). These results are largely confirmed by the analysis of the 
He 1 A 5876 line in the April 2005 SPM time series. Whilst the 
variations in the wings of the line are dominated by two low fre- 
quencies at 1.08 and 1.30d _1 , v\ (or its alias at 12.59d~') as well 
as V2 are detected with a roughly constant amplitude over the en- 
tire width of the line profile and the blue to red phase differences 
amount to ~ 7 n and ~ 4.5 n for v\ and V2 respectively. 

The SPM echelle spectra also cover the He 1 A 447 1 line. The 
frequency content of the power spectrum of the latter essentially 
confirms the results obtained for the Hei A 5876 line: a domi- 
nant peak at 0.08 d -1 (alias of the 1 .08 d -1 peak reported above), 
as well as two pulsation signals found at vj as well as V2 - 1. 
The main difference with the above results concerns the blue to 
red phase differences which, in the case of the Hei A 4471 line, 
amount to ~ 5.8 n and ~ 3 71 for vi and \>i — 1 respectively. 

A similar description applies also to the April 2007 peri- 
odogram of the He 1 A 5876 time series. Indeed, the highest peak 
is again found at low frequencies, 1 .06 d -1 (along with its aliases, 
Fig. [9]). Prewhitening the time series for this frequency yields 



secondary peaks at 8.25 (= V2) and 14.61 d -1 (likely the v\ + 1 
alias). Prewhitening the frequencies 1.06, 8.25 and 14.61 d , we 
find that the blue to red phase difference amounts to 6.5 n for 
the 14.61 d -1 frequency (the same result holds if we adopt v\ in- 
stead). For V2 we find a phase difference of about 4 ji (see Fig.[9j>. 

3.3. Photometry 

We have analysed the time series of HD 93521 in each of the 
seven filters of the Geneva system using the Fourier method of 
Heck et al. ( 119851 I. For each filter, we find that the periodogram 
indicates significant power at frequencies below about 10 d~'. 
In all cases, we find that the power spectrum can efficiently be 
prewhitened with three frequencies (see e.g. Fig.lTOt. 

Some frequencies are found consistently in several filters 
though with different amplitudes. This is the case for 7.90 d 1 
(found in six filters out of seven), 0.88 d -1 (detected in four fil- 
ters) and 2.00 d 1 (detected in three filters). None of the frequen- 
cies reported here is observed in the April 2005 spectroscopic 
time series that was obtained in coordination with the photomet- 
ric campaign. 

The U photometry yields somewhat different results from 
most other filters: the variations are clearly dominated by the 
2.00 d _I frequency, while in most other filters the various fre- 
quencies appear with similar amplitudes. At first sight, this 
might seem surprising because, for OB stars, the signatures of 
non-radial pulsations are usually expected to be quite strong in 
the U filter (De Cat et al. 120071 1. We caution that frequencies 
close to an integer number of d _I (such as the 2.00 d _I fre- 
quency) could actually be artefacts due to the interplay between 
the sampling of our time series and a few deviating data points. 
In this respect, we emphasize that the U filter is most heavily af- 
fected by atmospheric absorption and therefore, the 2.00 d -1 sig- 
nal could be related to imperfect extinction corrections, at least 
in the U filter. 
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February & April 2006: He I X 5876 
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Fig. 8. Left: same as Fig.[4]but for the He i A 5876 line as observed in February and April 2006. The power spectra are averaged over 
the wavelength domain from 5860 to 5900 A. The frequencies used for prewhitening are 1.10, 12.63 and 8.26 d -1 . Right: same as 
Fig.|6]but for the Hei yl 5876 line as observed in February and April 2006. The semiamplitudes and phases of the 1.10, 13.60 and 
8.26 d -1 frequencies are shown from top to bottom in the right panel. 
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9. Same as Fig.[8]but for the He i A 5876 line as observed in April 2007. The frequencies used for prewhitening are 1 .06, 14.62 
8.25 d" 1 . 
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Table 2. Results of the Fourier analysis of the photometric time series of HD 93521. Columns 2, 3, 7, 8, 9 and 10 yield the average 
magnitude and the dispersion about the mean for HD 93521 and the two comparison stars. Columns 4, 5 and 6 list the highest peaks 
in the periodogram of the photometric time series of HD 93521. The numbers in brackets provide the semiamplitude of the variation 
associated with each of these frequencies. 



Filter 








HD 93521 






HD 90250 


HD 96951 




mag 


cr 






Frequencies 






mag 


cr 


mag 


cr 


V 


7.0043 


0.0047 


7.90d- 


(2.9 mmag) 


2.05 d" 1 (2.5 mmag) 


7.10 d~ 


(2.3 mmag) 


6.4900 


0.0032 


7.8702 


0.0035 


B 


5.7589 


0.0053 


0.88 d" 


(3.4mmag) 


7.90 d _1 (3.2 mmag) 


6.11 d - 


(2.4 mmag) 


6.8883 


0.0045 


6.9591 


0.0038 


U 


5.8809 


0.0079 


2.00 d" 1 


(11.2 mmag) 


6-lOd" 1 (3.5 mmag) 


7.90d- 


(3.4 mmag) 


9.0681 


0.0085 


8.4939 


0.0060 


Si 


6.5113 


0.0057 


0.88d- 


(4.9 mmag) 


4.03 d" 1 (2.8 mmag) 


7.90 d" 


(2.6 mmag) 


8.2339 


0.0052 


7.8505 


0.0048 


B 2 


7.3870 


0.0057 


0.88 d~ 


(4.4 mmag) 


7.90 d -1 (4.3 mmag) 


6.81 d~ 


(3.1 mmag) 


8.0003 


0.0039 


8.4397 


0.0047 


V, 


7.6856 


0.0057 


2.00 d" 


(5.4 mmag) 


0.88 d 1 (4.3 mmag) 


7.90 d" 


(3.0 mmag) 


7.2655 


0.0040 


8.5661 


0.0050 


G 


8.2383 


0.0056 


0.04 d" 


(4.3 mmag) 


7.58 d" 1 (2.8 mmag) 


4.93 d" 


(2.5 mmag) 


7.4440 


0.0036 


9.0553 


0.0055 
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Fig. 10. Fourier power spectrum of the B\ data of HD 93521. 
The top panel yields the raw power spectrum, whilst the mid- 
dle panel corresponds to the power spectrum of the time se- 
ries prewhitened for three frequencies (2.04, 6.74 and 8.54 d _1 ). 
The lower panel displays the spectral window of the photometric 
time series. 



Table|2] also lists the mean magnitudes and the l-cr stan- 
dard deviations of the measurements of HD 93521 as well as of 
two reference stars (HD 90250, Kl III and HD 9695 1 , Al V) that 
were observed during the same campaign. The dispersion of the 
measurements of HD 93521 is always larger than those of the 
reference stars including in the filters where HD 93521 is signif- 
icantly brighter than the reference stars. The only exception is 
the U filter where the largest value of <x is found for HD 90250, 
however this star is by far the faintest in the U band. We fur- 
ther note that the Fourier power spectra of both reference stars 
are quite different from those of HD 93521 and are essentially 
consistent with white noise. 

In summary, HD 93521 appears to be slightly variable in 
photometry, though the frequency content of the variations is 
rather complex and the frequencies found in the photometric 
time series have very small amplitudes. It is likely that an even 
longer time series of very accurate photometric measurements 



(maybe from a space-borne observatory) is needed to disentan- 
gle the frequency spectrum. 

4. Discussion 

4. 1 . General properties 

A summary of the amplitudes of all the frequencies detected in 
the spectroscopic part of this study is provided in Table|3] 

A first conclusion from this table is that the v\ frequency is 
not always detected: while it clearly dominates in the April 2005 
data of the He i A 6678 line, it is absent from our February 1997 
time series of the same line. Also, this frequency is never de- 
tected in the Ho- line profile variations. In the Hei /U5876 and 
6678 lines, this frequency is associated with a maximum semi- 
amplitude of modulation of 0.004 - 0.006 in units of the con- 
tinuum. On the other hand, V2 or its aliases are detected in each 
line and each observing campaign. The maximum semiampli- 
tudes of this mode are 0.003 (Hei ,U5876, 6678 in 2006 and 
April 2005), 0.006 (Ha, February 1997) and 0.008 (Hei ,16678, 
February 1997) in units of the continuum. 

Another conclusion is the lack of any significant signal at 
the V3 = 2.66 d -1 frequency (P3 = 9.0+ 1.2hr) that was reported 
by Howarth et al. ( 1998 ) from their analysis of 103 IUE spectra 
with a median sampling of 0.60 hr. This frequency is thus appar- 
ently absent from our data, except perhaps for the detection of 
the 2.51 d" 1 signal in the April 2005 SPM Ha data. 

Generally speaking, the (low-level) periodic variations with 
frequencies v\ and V2 are not the dominant source of line profile 
variations in any of the lines investigated here (except for the He 1 
A 6678 line in April 2005). In the majority of the cases, the most 
important variations are actually characterized by low frequen- 
cies. These low-frequency variations affect both the core of the 
lines as well as the wings and emission lobes (in the case of the 
Hot and He 1 A 5876 lines). Apparently, these modulations do not 
occur with a single stable clock: none of the low frequencies is 
detected in more than one line and for more than one observing 
campaign (except perhaps for the frequency around l.lOd" 1 ). 
It seems thus more appropriate to talk about time scales than 
periods for these longer term variations. This result casts doubt 
on a rotational modulation as the origin of the long-term vari- 
ations: Howarth & Reid (1993) reported the analysis of 21 op- 
tical echelle spectra of HD 93521 acquired over two nights in 
February 1992. From the variations of the emission wings of the 
Hot and Hei A 5876 lines, these authors estimated a rotational 
period of 35 hr (v = 0.69 d -1 ). None of the lines investigated 
here displays a significant signal at this frequency. The closest 
detections are found for the 1997 OHP Ha data (the alias of the 
highest peak at 0.70d~') and for the same line observed in 2005 
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Table 3. Maximum semi-amplitudes (in units of the continuum) of the frequencies detected in the 2-D Fourier analyses of the line 
profile variability of HD 93521. For each line and each epoch, the amplitudes are given for those frequencies that were detected 
with a good significance above the noise level. 
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SPM04/05 


SPM04/05 


2006 
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~ 1 + V] 
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13.68 


V\ 
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0.006 
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(the alias of the highest peak at 0.78 d~'). Therefore, our data do 
not confirm the existence of a rotational period of about 35 hr 
and no unambiguous rotational period can be identified observa- 
tionally for this star. 

In the remaining subsections we will focus on the interpreta- 
tion of the variations seen at the v\ and V2 frequencies. 

4.2. vi and vt as non-radial pulsations 

In this section, we assume that the line profile variability at the 
frequencies v\ and V2 is due to two different non-radial pulsa- 
tion modes. For multi-mode pulsations, one expects to observe 
variability with the genuine pulsation frequencies and their har- 
monics, as well as with their sums and beat frequencies (see e.g. 
the c ase of the 09.5 V star £Oph, Kambe et al. [19971 Walker 
et al. 120051 1 . In HD 93521, we find no significant signal at the 
first harmonics of v\ and V2 nor at the sum or beat frequencies. 
Schrijvers & Telting ( 119991) argue that for non-radial pulsations 
with an amplitude of order 10% of the mean line depth, the ab- 
sence of a first harmonic is an indication that the line profile 
variability is mainly due to temperature effects (rather than to the 
Doppler-redistribution of flux). In HD 93521, the amplitudes of 
the vi and V2 modes in the He i A 6678 line are of order 5 - 10% 
of the maximum line strength. While these amplitudes might be 
somewhat low for the harmonics to be detected, we nevertheless 
note that temperature effects could play a significant role in the 
line profile variability observed in HD 93521. 

Telting & Schrijvers (fl997l l and Schrijvers & Telting (119991 
derived linear formulae relating the observable phase differences 
between the blue and red line wings to the degree I and the ab- 
solute value of the azimuthal order |m] of the pulsation mode. 
These relations are applicable to non-zonal (m + 0) spheroidal 
and toroidal modes and are valid also for multi-mode pulsations 
including those cases where temperature effects dominate over 
radial velocity effects. The blue-to-red phase difference of the 
mode yields /, whilst the phase difference for the first harmonic 
leads to the value of \m\ (Telting & Schriivers |1997l ). 

For the vi pulsation mode of HD 93521, the phase <f>\(A) is 
a monotonic function of wavelength (see Figs. [6] [8] and [9]), i.e. 



there are no huge changes in the slope over the interval where 
the significant line profile variations are detected. Therefore, this 
mode is unlikely to be an 'outlier' in the sense defined by Telting 
& Schrijvers (1997) and the relation between / and the blue-to- 
red phase difference A</> inferred by the latter authors should thus 
be applicable to these modefl Applying these relations to the 
phase differences given in TableH] yields / values of 8 + 1 for vi 
and 4+1 for V2. 

Howarth & Reid d 1993b interpreted the vi modulation as the 
signature of sectoral mode non-radial pulsations with I = —m - 
9. They further inferred a horizontal to radial velocity variation 
amplitude of k < 0.3. Later on, Howarth et al. (1998) derived 
/ =s 10+1 and 6 ± 1 for vi and V2 respectively, with m + I < 2. 
Whilst our values of I are in rough agreement with those of 
Howarth & Reid ( fl993l and Howarth et al. JT9981 , the lack of 
a significant power in the first harmonics prevents us from de- 
riving the value of \m\ from the simple scaling relations for this 
frequency. We note that the low amplitudes (or actually the upper 
limits on the amplitude) of the first harmonics compared to those 
of the genuine frequencies, suggest indeed that |m| is likely close 
to I (Schrijvers et al. l 1997b . However, we stress that in a rapidly 
rotating early-type star, such as HD 93521, the combined effects 
of a concentration of the pulsation towards the equator and of 
a non-uniform temperature distribution due to gravity darkening 
leads to more complicated amplitudes as a function of wave- 
length (e.g. Townsend 1997). A detailed line profile modelling 
is therefore required to derive the values of m and we defer this 
to a forthcoming paper. 

A priori, the shape of the semiamplitude of the modes as a 
function of wavelength (see Figs.|6j [8] and |9j) suggests that the 
modes have a rather modest ratio between the amplitudes of the 
horizontal and radial velocity variations {k < 0.5). However, 
Schrijvers & Telting (1999) showed that the typical double- 
peaked shape of the amplitude found for velocity-dominated line 
profile variability with high k values vanishes when temperature 



4 Telting & Schrijvers ( 1997 1 caution however that their relation was 
established for models with moderate rotation and hence stars that are 
not significantly flattened by rotation. 
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Table 4. Properties of the frequencies detected in the 2-D Fourier analyses of the line profile variability of HD 93521. The / values 
were derived from the blue-to-red phase differences discussed in Sect. |3.2| using the formula of Telting & Schrijvers ( 119971 ). 



y(d-') 


P(hr) 


Detection 


A(f>/n 


/ 


Comment 


14.61 


1.64 


He i ,15876 (2007) 


6.5 


7.2 


likely alias of V\ 


13.68 


1.75 


He i ,16678 (04/05) 


7.0 


7.7 


= v\ 


12.63 


1.90 


He i A 5876 (2006) 


7.0 


7.7 


likely alias of Vi 


8.31 


2.89 


He i A 6678, He i A 5876 (2006,2007) 


4,3.5 


4.5,3.9 


= v 2 


7.30 


3.29 


He i ,16678 (02/97) 


3.0 


3.4 


likely alias of v 2 



effects become important. Therefore, since we cannot exclude 
that the pulsation modes in HD 93521 might be affected by tem- 
perature effects (see above), we cannot make a secure statement 
about k. 
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Fig. 11. Top panel: normalised Fourier power spectrum of the 
He i A 6678 line as observed in April 2005. Middle panel: nor- 
malised Fourier power spectrum of the B\ photometric data ob- 
tained in April 2005. Bottom panel: product of the two nor- 
malised power spectra. 



Our photometric data do not reveal significant variability at 
the vi and V2 frequencies. However, because the v\ mode is of 
high degree, the integrated flux variability due to these pulsa- 
tions is indeed expected to be rather modest a priori and it could 
be that these variations just have too low an amplitude to be de- 
tected. To see whether or not this explanation is plausible, we 
have used the method outlined in Aerts et al. (2006). First, we 
have normalized the periodogram of the April 2005 He i A 6678 
time series by dividing it by the power of the highest peak. The 
same procedure was then applied to the periodograms (for each 
filter) of the photometric data, also taken in April 2005. The nor- 
malized periodograms of the spectroscopic and photometric data 
were finally multiplied two by two with the rationale that fre- 
quencies that are present in both data sets should be dominant 
in the product of the periodograms (see Aerts et al. 2006b . In 
this way, we find that the V2 frequency is clearly seen in all the 
products, whilst it had a significantly lower amplitude than v\ in 



the spectroscopic periodogram (see e.g. the case of the B\ fil- 
ter in Fig.rTTTi. The vi frequency is seen with a strength slightly 
larger than the V2 frequency only in the product of the spectro- 
scopic and photometric periodograms of the U and G filters. The 
fact that V2 appears more prominently in these products than v\ 
would be consistent with our conclusion that the former mode 
has a lower / value than the latter. Whilst this test is not a proof 
for the presence of the spectroscopic frequencies in the photo- 
metric variations, it nevertheless suggests that obtaining an ex- 
tensive photometric time series with a lower noise level (e.g. us- 
ing a space-borne observatory) would definitely be worth the ef- 
fort. 

We have further found several other possible periodicities 
that could be present in the photometric data. If real, these mod- 
ulations might correspond to either radial or lower degree non- 
radial pulsations. Such modes are difficult to detect in rotation- 
ally broadened line profiles, but might well produce an observ- 
able signature in the photometric data. Walker et al. (2003) re- 
ported on MOST photometry of ( Oph: they found that the light 
curve is dominated by a 4.633 hr period with a semiamplitude 
of 7.3 mmag, whilst the other modes have semiamplitudes be- 
low 2.25 mmag. These authors suggested that the light varia- 
tions of £Oph mainly result from radial pulsations. In the case 
of HD 93521, the vast majority of the frequencies detected in 
the periodogram of the time series have semiamplitudes that are 
significantly lower than what was found in the case of £Oph. 
Nevertheless, HD 93521 would be a very interesting target for 
an intensive photometric monitoring from space. 



4.3. The need for an alternative interpretation? 

Whilst multi-frequency non-radial pulsations offer an attrac- 
tive interpretation of the line profile variability observed in the 
spectrum of HD 93521, we must nevertheless ask the question 
whether there could be alternative explanations. We stress that 
the stability of the frequencies over many years (when they are 
detected) likely implies that the profile variations are produced 
by one or several stable clocks such as pulsations (considered 
above), rotation, orbital motion... Generally speaking, the dis- 
tinction between these different mechanisms is a non-trivial is- 
sue (see e.g. the discussion in Uytterhoeven et al. l2005b . For in- 
stance, in the case of ( Oph, Harmanec (1989) proposed a dif- 
ferent scenario where the moving bumps in the absorption lines 
actually stem from rotating inhomogeneities of the circumstellar 
material rather than from non-radial pulsations. 

In our case, there are also a number of reasons to consider 
alternative scenarios. Indeed, it is clear that currently a theo- 
retical model to predict the line profile variations produced by 
non-radial pulsations in a rapidly rotating massive star such as 
HD 93521 is still lacking. Therefore, the interpretation of these 
features within the framework of the available models requires 
an extrapolation that might be difficult to justify a priori. 
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Another feature that is certainly puzzling is the fact that the 
line profile variability is significantly detected only in lines that 
are potentially affected by residual emission possibly associated 
with a circumstellar disk (or flattened wind). Whilst the lack of 
a TVS signal in the purely photospheric O in A 5592 line (see 
Fig-UJ might be interpreted as this line forming near the hotter 
poles of the star where the pulsational amplitude is lower (as 
for He ii lines), the same explanation cannot hold for the photo- 
spheric features that occur at temperatures and gravities that are 
typical of those of the equatorial region (see the C n, N n and Si in 
lines discussed in Sect. 13. U . The latter features produce at best 
a marginal signal in the TVS (see Fig.[T}. However, regarding 
the possibility that the variability stems from rotating features 
in the circumstellar material, we first note that the frequencies 
vi and V2 are not detected in the emission wings of the lines 
analysed in this work. If the line profile variations were coming 
from the flattened wind, one would also expect them to affect 
the emission wings. We further note that all the lines where we 
have detected line profile variability are quite strong. Actually, 
the rather low amplitude of the profile variations might render 
them undetectable in the weaker equatorial lines. 

In relation to this, we note that a recent study of HD 60848 
revealed evidence for the existence of rather short (3.51 and 
3.74 hr) periods in the radial velocities derived from the emission 
lines of this 09.5 IVe star (Boyajian et al. 2007). These authors 
argued that these features might result from changes in the disk 
density or illumination caused by non-radial pulsations in the 
underlying star. Therefore, it seems that a disk origin for short 
period variations cannot be ruled out a priori, although it would 
also reflect the existence of pulsations in this specific example. 

Since HD 93521 was considered as a possible runaway ob- 
ject (Gies 1987J but see also the discussion in Sect.Q]), it could 
host a compact companion which could then trigger a periodic 
structure in the inner part of the disk that is seen projected 
against the photosphere. However, the presence of a compact 
companion should make HD 93521 a bright X-ray source, at 
least episodically when the compact object crosses the equa- 
torial wind (similar to Be-type high-mass X-ray binaries). We 
thus checked the ROSAT All Sky Survey images as well as sev- 
eral other X-ray catalogues. To the best of our current knowl- 
edge, there is no indication that HD 93521 is, or has ever been, a 
bright or even moderate X-ray emitter. Another point is that the 
detected frequencies v\ and V2 are much too high to correspond 
to the orbital period (> 13.5 hr and likely of order several days) 
of such a putative companion. A compact companion scenario 
appears therefore rather unlikely. 

Another point concerns the fact that both frequencies are too 
high to correspond to the rotational period of the star. In fact, 
whilst many of the physical parameters of HD 93521 remain un- 
known or poorly determined, we can obtain a rough estimate of 
the rotational frequency by considering typical parameters of a 
late O-type star. Let us assume that HD 93521 is an 08.5 stafl 
with a polar radius of 8 R Q and a mass of 20 M Q . If the star is ac- 
tually seen equator-on, the resulting rotational frequency would 
be 0.81 d -1 . In general, we can say that for any reasonable as- 
sumption on the stellar parameters, we find that v lot < 1 d" 1 . 
Now, if the line profile variations actually stem from a regular 
pattern of moving spokes in the circumstellar disk, v\ and V2 
would not necessarily have to correspond to v lot , but could rather 
be some harmonics of the latter (see Uytterhoeven et al. 2005). 



5 The observational classification as an 09.5 star is likely biased to- 
wards later spectral types as a result of gravity darkening due to the high 
rotational velocity and the nearly equator-on orientation of the star. 



We would thus have to look for a "super" frequency such that 
vi = n\ v sup and V2 = n2V sup , where n\ and «2 would be integer 
numbers. Within the uncertainties of our period determinations, 
a candidate for such a frequency would be v sup = 2.75 d _1 (with 
«i = 5 and «2 = 3). Again, the latter frequency is much too high 
to correspond to the rotational frequency. The super-frequency 
could be in agreement with the likely value of the rotational fre- 
quency if n\ and «2 were muliplied by 3. However this would 
imply a large number (15 and 9) of co-rotating structures around 
the star which seems rather difficult to explain. 

In summary, we conclude that, all the alternative scenarios 
envisaged here fail in explaining the modulations at the v\ and 
V2 frequencies. Hence, despite some difficulties, the multi-period 
non-radial pulsations model remains currently the most plausible 
explanation for the line profile variations seen in HD 93521. 
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